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Abstract: We propose the operatorial Baxter's TQ-relations in a general form 
of the operatorial Backlund flow describing the nesting process for the inho- 
mogeneous rational gl(K\M) quantum (super)spin chains with twisted periodic 
boundary conditions. The full set of Q-operators and T-operators on all lev- 
els of nesting is explicitly defined. The results are based on a generalization of 
the identities among the group characters and their group co-derivatives with 
respect to the twist matrix, found by one of the authors and P.VieirapQ. Our 
formalism, based on this new "master" identity, allows a systematic and rather 
straightforward derivation of the whole set of nested Bethe ansatz equations for 
the spectrum of quantum integrable spin chains, starting from the R-matrix. 
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1. Introduction 

It has been noticed long ago that the mathematical structures behind the quan- 
tum integrable spin chains have many similarities with the theory of classical 
integrable systems, such as KP or KdV hierarchies. It goes of course not only 
about the obvious correspondence between the quantum integrable systems and 
their classical limits when, for example, the quantum transfer matrix of a quan- 
tum 1+1-dimensional system becomes the classical monodromy matrix of the 
corresponding classical Lax connection. There is a more striking "classical" fea- 
ture of the quantum integrability: The quantum transfer matrix represents a 
natural (spectral parameter dependent) generalization of the Schur character 
of a classical algebra [5] given by the so called Bazhanov-Reshetikhin (BR) 
determinant formula and, as such, it satisfies a certain Hirota bilinear finite 
difference equation, which appears in the quantum context as a certain fusion 
relation among the composite quantum states appearing in quantum spin chains 
as certain bound states ("strings") of Bethe roots 0H1I3I3E] ■ Similar, though a 
more complicated realization of Hirota discrete "classical" integrable dynamics 
has been observed in the context of the quantum (l+l)-dimensional QFT's, or 
sigma-models [6 ,81(1], an observation which appeared to be at the heart of an im- 
portant advance in the study of the spectrum of the AdS5/CFT4 correspondence 

mam. 

Hirota equation immediately brings us to the idea that quantum integrabil- 
ity, at least for certain quantities, such as transfer-matrix eigenvalues, can be 
viewed as a specific case of classical integrability and of the theory of classical 
tau- functions. Indeed, the character of a classical group, say gl{K), is noth- 
ing but a tau-function of the KdV hierarchy. It was proposed in [T] to view 
the quantum transfer matrix of a rational quantum Heisenberg-type gl{K\M) 
(super)-spin chain with twisted boundary conditions as a quantum, operatorial 
generalization of the character and to construct the transfer matrix (T-operator) 
by acting on the character in a given irrep, as a function of the group element 
(twist), by special group derivatives, called the co-derivatives. The formalism of 
co-derivatives has led to a direct proof of the BR formula [T] (see also [T^] ), and 
the basic underlying identity for the characters found in [I] seems to be just a 
new form of the KdV Hirota identity (the fact yet to be understood). 

In the present paper, we want to move even further in this classical inter- 
pretation of the quantum integrability and to generalize the basic identity of [1] 
to include the Baxter's TQ-relations into our formalism. This implies a natural 
definition of all Baxter's Q-operators, rather different from the one known in 
the literature [TIfl[(M ^ IT ^ IT ^frnfT%l fT!T| ■ and more generally, of the T-operators 
on all levels of the nesting procedure. This nesting takes a form of a Backhand 
flow, directly for the T- and Q-operators. Due to the fact that all of them be- 
long to a commuting family of operators, all these relations can be immediately 
transformed into the well known functional form, for their eigenvalues |20[|211 

1221231. 

Our main identity given in the next section offers an interesting alternative 
and a concise approach to the quantum integrability uncovering the whole struc- 
ture of the nested Bethe ansatz, from the R-matrix and the Yang-Baxter relations 
all the way to the nested Bethe ansatz equations, in the general operatorial form 
for all the intermediate quantities. 
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2. Transfer-matrix, co-derivative and TQ-relations 

We recall that the main object of our study is the transfer matrix of an inho- 
mogeneous quantum spin chairQ 

TW(«) = tr A (R$\u-G N )®---®R { l X} (u-e 1 )ir x (g)) (1) 

where ir\(g) is a matrix element of a twist matrix g £ GL{K) in an irrep A and 

R| A} (u) = uI + 2^ ef a <g> 7T A (e a/9 ) = u + 2 7\ A (2) 

a/3 

is the R-matrix in irrep A in auxiliary space (and in fundamental irrep in the 
quantum space). The identity operator I is implicit in the r.h.s., it will often 

be omitted for the brevity. The gl(K) generator corresponds to the i th 

quantum spac^l (which is in the fundamental representation) and ir\{e a p) - 
to the auxiliary space. When A is also fundamental then 7\ A becomes a usual 
permutation operator ef a = I®^ 1 ) Oe^igil®^-^ so that V t = J2 a p 4^®e Q(3 . 
Vi permutes the indices of the auxiliary space and the quantum subspace. This 
(u) is a spectral parameter «£C dependent operator on the quantum space 
W. It is polynomial in u and in the inhomogeneities ^ gC. 

The main goal is to find all the eigenvalues of this transfer matrix for the gl{K) 
quantum spins. For that we work out an operatorial Backlund transformation, 
which can be also called the nesting procedure, whose main goal is to derive, in 
a deductive way, without any assumptions, the nested system of Bcthc ansatz 
equations defining these eigenvalues. On the way, we will encounter a collection of 
the intermediate T-operators, and the Baxter's Q-operators as their particular 
case, at each level of nesting. The operatorial TQ-relations, representing the 
Backlund transformation reducing the problem for a gl(k) subalgebra to a similar 
problem for the gl(k — 1) subalgebra, in the nesting procedure corresponding to 
the chain of embeddings gl(K) D gl{K — 1) D ••• D gl(l) were given in their 
functional form in [21] (for the super-spin chains in [20]). All the T- and Q- 
operators in this nesting procedure are labelled by the subsets I of the full 
set 3 = {1,2,..., K} as Tj A ^(u) and Q/(u). The original transfer matrix ( |2"]l 

corresponds to the T-operator for the full set (it) = T^(u). There are 2 K 
subsets of the full set 3 and they can be described in terms of the Hasse diagram 
based on the inclusion relations. A chain of subsets of the full set 3 = Ik D 
Ik-i D • • • D 7q = 0, where |Jfc| = Card(Zfc) = k, forms a path on the Hasse 
diagram. We will call this the nesting path, and (K — k) - the level of nesting. 
There are K\ different nesting paths for gl(K) (see fig llb[) . This description of 

1 Throughout this paper, all operators will act on the same Hilbert "quantum" space 
H = (C K )® JV (resp (C^ ' M )® iv for supersymmetric spin chains). We use the Young diagram 
A to label the irreducible tensor representation. The fundamental representation corresponds 
to the Young diagram with one box. 

2 Here {e,j}|^_ 1 satisfy the commutation relations [ejj,e&;] = S-jk e il ~ ^li e kj- The i th 

quantum space is the i th factor in H = (C K ) <X> (C^ ) (g) • • • ® (C^). In this way. e.g. ei^ = 

v ^ , 

N times 

1) (gi e ^ a (g,J®( N — l ). The generalization of our construction to the case of gl(K\M) super- 
spins will be given in secl5l 
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all the 2 K Q-operators based on the Hasse diagram was proposed in [23], and 
will be used throughout this paper. 

The general TQ-relations, derived in sec|4]in the operatorial form, are given 

by 

Tf («)Q/j(«) = TJ ,jHQ/(u) - XjT^iu + 2)Q 7 ( M - 2), (3) 

where g — diag(xi, x%, ■ ■ ■ ,xk) is the twist matrix in the diagonal basis, the 
superscript s in T-operator denotes the symmetric A = (s) irrep in the auxiliary 
space, by I C 3 = {1, 2, • • • ,K} we denote a subset of the full set of indices 
(labeling the eigenvalues) and by I,j = I U {j} C 3 we denote a subset with 
one more index j ^ /. This TQ-relation relates the T-operator Tf and the T- 
operator Tf ■ of the previous level of nesting (which has one more index) . A 
chain of these relations allows to relate the original transfer matrix TJ(tt) on 
the level zero of nesting to the u- independent operator T2(m) given by (|4"Tj) . We 
will also show that the Q-operators are equal to the T-operators taken at an 
empty Young diagram: 

Q/(u) = T?(«). (4) 

In the papers [2"T1|20| ■ all the T- and Q-operators at intermediate steps were 
assumed, by self-consistency and without a proof, to be polynomials in u. This 
analyticity assumption immediately leads to the nested Bethe ansatz equations 
defined by the nesting path. In this paper, we complete the missing link of the 
chain and find the explicit operatorial form of the Backlund flow ([3|). 

In what follows, we will extensively use the definitions and the identities of 
[TJ. In particular, the co-derivative D defined there and used through the whole 
current paper is a very simple object defined by its action on any function of g 
as follows 

D®f{g) = -^-®f{^ e g) (5) 

d <t> 0=0 

where is a matrix in the fundamental representation: <f> ■ e = Y] a a Sap^ an( i 
^ = J2 a /3 ea P~Q^F- ^ s mam property, which also could serve as its definition, 
manifests in its action on the group element in fundamental irrep: 

D®g = V{l®g) 

where V is the operator of permutation between the 1 st and the 2 nd spacefl 
Many other useful properties of this co-derivative, mostly following from the 

application of the standard Leibniz rule can be found in pQ and some of them 

are summarized in the appendix IA.1I 

Using the co-derivative we can for example rewrite the T-operator (|TJ) in the 

following way: 

T {A} (u) = (ui + 2D) <g> (ua + 2D) <g> • • • ® (u N + 2D) X{A} (g) (6) 



3 Explicitly in indices, the last relation looks like D'.J g % ? = . It is a usual matrix 

derivative obeying the Leibniz rules. Throughout this paper, we define the tensor (or matrix) 
indices of any operator A with respect to a basis {v^ l ®v k 2 <£>■ ■ ■CS)i'fcj V } as Aw^ ®vi 2 ®- ■ ■®'Vl N = 

Efej k 2 ■■■ *;« ^h 'l2*'-'" iff v ki ® v k 2 ® • • ■ ® v k N ■ 111 particular, a usual notation of a matrix 
element A^ i is written as A*. 
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(a) Hassc diagram for the Q-operators: gl(3) 
case. There are 2 3 = 8 Q-opcrators. 



Q1234 




Q123 Q124 Q134 Q234 




Ql2 Ql3 Q14 Q23 Q24 Q34 




Qi Q2 Qs Q4 




Q 



(b) Hasse diagram for the Q-operators: gl(4) case. There are 
2 4 = 16 Q-operators. 

Fig. 1: Hasse diagram for the Q-operators in gl(3) and gl(4) case. In (b), the 
red thick lines denote one of the 4! = 24 nesting paths characterized by a chain 
of the index sets 3 — I4 D I3 D 1% D I\ D Iq, where I4 = {l,2,3,4},/ 3 = 
{1, 2, 3}, I2 — {2, 3}, I\ = {2}, I — 0. Note that (b) contains the diagram (a) for 
gl{2>) as a subdiagram. 
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where itj = u — Oi and x {A} (3) = tr 7T\(g) is the character of the twist g in the 
irrcp A. The action of each of N brackets adds a new spin to the system, with 
its fundamental quantum space. 



3. The master identity 



In this section we will formulate the main identity of this paper, called here the 
master identity (|8|) - the basis of our approach to the quantum integrability. It 
is proven in the appendix [El 

This identity involves the generating function of characters in symmetric (A = 
(s)) or antisymmetric (A = (l a )) tensor irreps, 



w(z) = det 



1 



1 



1 ~ z 9 UU0--ZXJ) 



1 



Er= (-i)^ a x (a) (ffr 

where the eigenvalueaj of the twist matrix g, and z G C. This master 

identity will relate operators of the form <E>i(ui + 2D) n,k, u '(^fc) (^ or an arbitrary 
set of complex numbers {£&}), which is a generalization^ of the T-operator ©. 
To avoid the bulky notations in this definition, we assume all the terms like Ui 
and 2+Ui to be multiplied by the identity operator I, and the tensor product (g>i is 
taken at0 unless it is explicitly stated otherwise. Due to the commutation 

relation (I7ip following from the Yang-Baxter relation, the operators <S>i(ui + 
2D) Y\k w (tk) are conserved charges, in the sense that they belong to the same 
family of commuting operators as the T-operators. 
The master identity reads as follows (t, z G C): 



(t-z) (E> l (2 + u l + 2D) w(z)w(t)n ■ ®i(ui + 2D)II 



= t <g>i(ui + 2D) w(z)II ■ <E>i(2 + iii + 2D) w{t)II 



-z <g>i(2 + m + 2D) w(z)II ■ ®i(v,i + 2D) w{t)II , (8) 

and it holds for any function 11(g) of the form 11(g) = Y[k w (^k), or equiva- 
lently, 11(g) = det(f(gj) — J\j=i f( x j) where f(z) is an arbitrary fixed function 
analytic in the vicinity of z — 0. In (JSJ), the dots between consecutive brackets 
stand for multiplication of operators acting on the quantum space, and each 
co-derivative operator D acts on what lies to its right inside the square brackets. 

The proof of our main identity © is given in Appendix [El but it can be easily 
proved directly, for a few small N's, on Mathematica. The identity represents a 
natural generalizatiorQ of the equation (4.1) [equation (20) in the arXiv version] 

in DP- 



4 Throughout this paper, we shall assume that Vi 5^ j , Xi 7^ xj. 

5 The T-operator in symmetric irrep can be indeed obtained as T s (ti) 
^(ll) S [®l(«i + 2D) Mz)]^ Q 

6 ~SfL 1 Ai = Ai Cg) A2 ® • • • SS Ajv for any indexed operators {Ai}^_ 1 . 

7 The identity (4.1) of QQ corresponds to a particular case when 77 = 1 and Vi, Ui = 0. 
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To conclude this section, let us demonstrate the use of the master identity 
([8|) by deducing from it a particular case of the Hirota relations for the transfer 
matrices in particular representations. For that purpose, let us focus on the case 
when II = 1. Then, by expending w's in symmetric characters according to (J7J 
and keeping the coefficient of t s z s in each term of (JEJ), we get the relation 



,•(2 + m + 2D) (xs-iXs' - XsXs'-i) 



(Ui + 2D) X s> ■ ®i(2 + Ui + 2D) xs-i 



® i (2 + w J + 2L>) X s'-1 



i(Ui + 2D) X s 



(9) 



Choosing s' = s + 1, and rewriting (jH]) usingjf] ^ , we get a relation in terms of 
the T-operators of rectangular Young diagrams : 

-T( 2 < s >(m + 2)T<°''>(u) = T^' s+1 \u)T^ s -^(u + 2) - T^' s \u + 2)T< 1 '">(u) 

(10) 

where (a, s) denotes the representation with an a x s rectangular Young diagram 
A = (s a ), i.e. T^ s \u) = T« sa »(M). 

Eq. (fTOf , is a particular case of general Hirota relations (|35|) for the fusion in 
rational spin chains known since long and proven in [I] at zero level of nesting 
(1 = 3). 



4. Baxter relations for T- and Q-operators 

In this section we will derive from our main identity the operatorial Backlund 
flow in the form of the TQ-relations described above, and even more generally, of 
TT-relations at every step of the nesting procedure, as well as the QQ-relations 
|24|I16U25P 17 P20U22P26U18] (see also earlier papers [27]. and a recent presentation 
in [53] used in this paper) which give an immediate access to the full set of nested 
Bcthe ansatz equations (also written in an operatorial form in quantum space in 
the subsection l4.5p . At the same time, it will give a natural operatorial definition 
of these quantities on every step of the nesting, and in particular of all the 2 K 
Q-operators. Since all these T- and Q-operators belong to the same family of 
mutually commuting operators, we can transform these relations, at any stage 
of the nesting procedure, to the operatorial ones, for T- and Q-operators. 



4-.1. First level of nesting. Now we will obtain from the master identity ([5J the 
operatorial Baxter's TQ-relations. We will start from the first level of nesting. In 
what follows we will frequently use the notation / = 3\I for the complimentary 
set of /. In particular for any element j € 3, we use a notation j = 3 \ {j}. 

8 We also use the fact that the characters for rectangular representations y( a ' s ) = ^ , 

satisfy a simple Hirota relation x (o ' s+1) x (a ' s_1) - x (a,s) X (a,s) = -x (a_1 ' s) x' a+1,s) ■ 
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Definition of Q- operators. 

In accord with (HJ and © (where T^ A ^(u) denotes (u)), the Q-operator 
on the zero level of nesting is, by definition, 



Qx2...k(u) = Qfl(w) 



' N > > 



|«>JV 



(11) 



which is a simple function of Ui's, times the identity operator in the full quantum 
space. In particular, the last factor in the l.h.s. of (JSJ) becomes Q^w) when 
77 = 1. 

We will see in what follows that the Q-operators of the first level of nesting 
Qj(u) can be defined through the residues at the poles in the expression: 



(i-gtf 



) N 



)( Ul + 2D + 2) w(t) 



= } — h polynomial in t, (12) 

Z ' 1 — Xnt 



where the normalization factor (1 — gif® 1 ^ is necessary in order to have only 
simple poles. Indeed, the co-derivative acting on w(t), having simple poles at each 
t = xj , j = 1, 2, . . . , K, produces the double poles in the same points, as it is 

clearly seen from (|64| in the diagonal basis. The factor (1 — <?t)® ^ transforms 
them again into simple poles, thus justifying the pole expansion in the r.h.s. of 

mi. 

The equivalent definition of the Q-operators is 



Qj(u) = lim (1 - xjt) (1 - gtf 



)N 



)(2 + Ui + 2D) w(t) 



(13) 



This Q-operator acts on the same quantum space H = (£, K )® N as the T-operator 
(J6j) . It is also important to notice that the Q-operator Qj(u) looses its dependence 
on some Ui's (see appendix |B1). For instance, for one spin (N — 1), if we denote! 
by |e&) a basis of eigenstates of g, we get 



Qj(w)|e fe ) 



((ui + 2) (1 - Xk/xj) + 2x k /xj) ILej i-xt/x, \ e k) tfk=£j 

(14) 



where we see that on the state | ) , Qj(w) is it-independent, while the action of 
other states is less trivial. 



T-operators and TQ-relations. 

Now we will transform the master identity ([8]) into a set of TQ-relations J3J 
on the first level of nesting]. For that we simply put II = 1. Multiplying © by 

9 In this definition of the "diagonal basis", | } is such that g\ei c ) = Xk\ef z ). 
10 The "first level of nesting" means that we will relate the original T- and Q-operators 
labeled by the full set = C! with some T and Q-operators labeled by j, which has one index 
less. 
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the matrix (1 — gt)^ N which commutes with all the factors of both the L.H.S. 
and the R.H.Slij, and picking the poles at t = l/xj we come to the equation 



lim (1 - Xjt)(l - z/t) (1 - gif 



®i(2 + m + 2D) w(z)w(t) 



Qi(«) 



ii(ui + 2D) w(z) •Qj(u)-2Cj ®i(2 + Ui + 2D) z w(z) -Qj(u-2). 

(15) 



It is useful to note that the factor (1 — z/t) ~ (1 



in the L.H.S. can be 



carried over to the right of the co-derivatives D allowing to use the relation 



(1 — z Xj)w(z) = (1 — z Xj) det — 



1 



det- 



1 



z 9 



z 3t 



where g- — diag(a;i , x%, ■ ■ ■ , Xj-i, Xj + \ , • • • , xk) in the diagonal basis. The possi- 
bility to move this factor across the derivatives comes from the factor (1 — gt)® N 1 
introduced to avoid poles of higher orders in (|12j) . Indeed, for example in the 
simplest, one spin case N = 1, we can easily check that (xjl — g) ■ D Xj = 
(xjl — g)Pj ■ Xj — x 2 jPj — gPjXj — 0, where Pj is the projector on the j th 
eigenspace of g. The generalization to any iV is rather trivial and is discussed in 
the Appendix IFl 

Now we introduce the characters of the first level of nesting 
Xs(ffj) = X{\=(s)}(dj) corresponding to the symmetric tensor representations 
of the sub-algebra gl(K — 1) c gl(K), defined by the generating function 



wj(z) = det 



1 - z a- rifcgji 1 - z x k) 



^2 zS Xs{gj) 



(16) 



and define the T-operators of the first level of nesting labeled by a Young diagram 
X=(s): 



Tf (u) = lim (1 - Xj t)(l-g tf 



)( Ul + 2 + 2D) Xs(9j)w(t) 



(17) 



The last formula also allows for an alternative to definition (|12l) of the Q- 
operators as of the T-operators for the singlet irrep in the auxiliary space corre- 
sponding to an empty Young diagram s = 0: Qj(u) = Ti =0 (u) confirming the 

(more general) relation announced in 

In the definition (|17j) , we take a residue of a given pole and at the same time 
use a character Xs (<?-) where one eigenvalue is "removed" . This "removal" will 
be at the heart of our nesting procedure, and its repetition defines a certain 
Backlund flow (nesting path). 



11 Which is clear in the diagonal basis since (1 — gt)^ 1 obviously commutes with permu- 
tations, and with tensorial product of diagonal matrices, and hence with any operator of the 
form ®i(a.i + D) w(b), due to its diagrammatic expansion given in appendix lA.il 
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Let us take the coefficient of z s (s e Z>o) in (fT5j) . analogously to what was 
done in ([9]), and rewritcF^l it, using (fT7]h as follows: 

Tf(u)Qg(u) = Tf(u)Qj( U ) - ajj-Tj- 1 ^ + 2)Q 7 ( U - 2), (18) 

which is the simplest Baxter's TQ-rclation in the operatorial form, the first of the 
chain of Backlund transformations among the commuting T- and Q-operators 
of the zero th and first level of nesting known for a long time [H] on the level of 
their eigenvalues. Here T|-(ii) = T^ s '}(it) is the T-operator of the zero th level 
of nesting, the original transfer matrix (JlJ, or (j6|), in the symmetric tensor irrep 
A = (s). The T- and Q-operators labeled by j = 3 \ {j} have K — 1 indices, and 
are considered in the first level of nestings- 
Let us also note that the T-operators can be also defined as the residues at 
the poles: 

(g)K + 2 + 2D) (t X {s}{g) - X{ s -i}(g)) w(t) 

i 

+ polynomial in t. (19) 

It is clear from these definitions and from (fTTj) that all these T- (and hence the Q- 
)operators, belong to the same family of commuting operators [T|(ii), T|, (u 1 )] = 

[Tj(u), Tg (it')] = 0. It will be also shown for all T- and Q-operators, on all levels 
of nesting. 



A t Ti (it) 



^ 1 

7 = 1 



Xjt 



)N 



4-2. Next levels of nesting . Now we will generalize this procedure, and the cor- 
responding TQ-relations, to all nesting levels. Suppose we want to consecutively 
"remove" the eigenvalues Xj t , Xj 2 , . . . , Xj k from the characters in the definition 
of T-operators, where / = {ji,j2, ■ ■ ■ >jk} is a subset of the full set of indices: 
I C 3 (their order is not important but they are all different). At such arbitrary 
level of nesting, we define a normalization operator 

B T =n(l-»i* J 0-(l-ff*i) 8JV (20) 
and the following product of generating functions of characters 

ii 7 =n w fe). (21) 

The definition of the Q-operator labeled by a subselE3 / = 3 \ I of the full set 3 
becomes 

where |I| = Card(T) =K-|I| 

(22) 

As explained above, we also use the fact that w(z)(l — xjz) = wj(z). 

In the same spirit, the fc th level of nesting will involve the quantities with K — k indices, 
the subset / defines the node on the Hasse diagram where the nesting process has arrived. 



Q/(lt) = lim By 



jei 



2D) n-, 
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and once again, it is an operator on the quantum space (C K )® N , which is poly- 
nomial (of degree N if I ^ 0) in the spectral parameter u. Its eigenvalues have 
degree < N, and it can be shown (see appendix |B|) that Qi(u)\ek 1: k 2 ,---k N ) is 
independent of all u n such that k n E I. 

We will show how to write a TQ relation between the T- and Q-operators 
labeled by I and the operators of the previous level of nesting, labeled by iUjk- 

Let us first generalize (|16p to define the characters of g I = diag ((xj)j e j): 



wj(z) = dct 



1 TT 1 a i \ w ( z ) 



(23) 



If we chose in the master identity 



77 = n 7UJZ = n 7\i k > 1 = and u ^ 
151 



u + 2| J U j k = U + 2\I\ - 2, then after multiplying it by 



t B 7 B, 



and 



taking the limit tj 
= lim 



we get 



5 i (2|7|+ Ui + 2£>) u;(*M*iJtfjqfc 
Di(2|J| - 2 + + 2Z>) w(z)n 7 



■ Qiuj k (u) 



'Uj fc 



Qi(u) 



B 



-00 



®i(2|I| + + 2.D) zwiz)!!^ ■ Qj{u - 2) , (24) 



where I U jk — {ii,j2, ■ • ■ ijfc-i}- These expressions are obtained by rewriting 
the z- independent factors using the formula ([22)1 . For instance, the last factor of 
the last term obtained from (151) is 



B 7 



®i(2|/| 



M + 2D) w(t jk )n Ah 



3,(21/1 



2D) n T , 

(25) 



which becomes Q/(lt — 2) when the limit tj — >■ — is taken. 

We define the T-operators for symmetric tensor representations as follows 



Tj(u) 



lim B 7 



N 



)( Ul + 2D + 2\I\) X s(9i)n T 



(26) 



where Xs(9i) is defined by (|23[) . Then the action of co-derivatives on Xs(<?j) is 
a priori rather nontrivial. The recipe to avoid this complication and compute 



15 Once again, the normalization factors By and B jrjj- commute with all the other factors, 

because they commute with all permutations, and with all operators g 11 ®g 12 [§>••■ g 1N , which 
are the building blocks of all other factors. 

On the other hand, we will see a posteriori in appendix El that the factor wj(z) can be freely 
move across the _D's. 



12 



Vladimir Kazakov, Sebastien Leurent, and Zengo Tsuboi 



T-operators is given in appendix iFl In terms of the generating function of Tf (u) 
we have 



Wi{u,z) = Vz s T*(u) = lim B 7 



s=0 



N 



,i=l 



( Ul + 2£) + 2|/|) ^-i7 7 

w 7 (z) 



(27) 



— 3-— lim By 

j&7 



N 



)( Ui + 2D + 2\I\) w{z)n T 



i=l 



Finally, we can notice that the L.H.S. of 
Then, expanding 



contains 



(28) 

(l-zx jk ) _ 1 



with respect to z and taking the coefficients of z s , we 
obtain using the definition (f2"5)) the following operatorial TQ-relation 

Tf( U )Q Juife (n) =TJ uifc (u)Q / («) - ^T|-)> + 2)Q z (u - 2). (29) 

It generalizes the similar, obvious relation among the characters of symmetric 
tensor irreps: 

Xs{gi) = Xs(gi,j) - XjXs-i(gi,j), (30) 

where j el and s G Z>i. 

Notice that again the Q-operator on any level of nesting is equal to the T- 
operator, with the same index set /, for an empty Young diagram A = (which 
corresponds to s = case): 



Q/(«) = Tf\u) 



(31) 



4-3. Generalization to any representations. There is a natural way to generalize 
the T-operators to any irreps A in the auxiliary space: 



t| A} (u) = lim By 



N 



)(u l + 2D + 2\I\)xMn T 



(32) 



where the gl(K — \I\) characters of the irreps A are given through the characters 
of the symmetric tensor representations Xs(9i) by the Jacobi-Trudi determinant 
formula 



X{X}{9i) = det Xx 1+t - 3 (9i) » 



l<z,j<a 



(33) 



where a is the number of rows in the Young diagram A. It is noteworthy that, due 



to the definition 



if A has more than |/| rows, because \{\} idi) = 0. 
The Bazhanov-Reshctikhin (BR) formula proven in [I] at the zero th level of 
nesting is also true for the T-operators on every level of nesting, and it reads 



-,{(Ai,A 2 ,-.Ao)} 



U.lZlQi(u-2k) 



det (T^ +l -\u + 2-2i)\. (34) 
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In the particular case of the rectangular Young diagrams A = (s a ) this means 
that the T-operators satisfy the immediate equivalent of the BR formula - the 
Hirota equation^, which is the same on any level of nesting, namely 

Tf ' s) (w+ l)T ( f' s \u- 1) = T^ 0+1 '* ) (u + 1)tS.°- 1 '* ) («-1) 

+ if ' s+1) [u - l)T^' s - l \u + 1). (35) 

It can be proven that it is a consequence of the equation (|86|P^I of appendix 
lEl It is also well known that the TQ-relation ([25)) implies the generating se- 
ries expression ((75)) of T-operators for symmetric representations in terms of 
Q-operators, detailed in the appendix [C] (see for instance [5T] at the level of 
eigenvalues). Then the Bazhanov-Reshetikhin determinant formula (or, equiva- 

lently, the Hirota equation) allows to write arbitrary t| A ^ (u) operators in terms 
of Q-operators and to check the following bilinear relations on the nested T- 
operators [H] (see |H] for the case with non-zero twist g): 

T (a+M) (u)T (a,s) (u) _ T (^)( U ) T ("+M) (u) 

- x J Tf+ 1 ' s ~ 1) (u + 2)T ( f' s+1} (u - 2), (37) 



T Ks+l) (u)T (a,s) (u) _ T M( U ) T (<M+1)( U ) 

- XjT^+^iu + 2)T ( I a - 1 ' s+1) (u - 2). (38) 

The TQ-relation ([29)) is a particular case of equation (|37|) when a = 0, and the 
two equations (|3"T1 |3"51) coincide with the definition of Backlund flow given in [2"2"j 
up to the permutations on the index set /, so that the definitions (|31[ I32p for 
nested T- and Q-operators explicitly give a solution of the linear system (|371 

EH). 



4-4- QQ-relations. In the previous subsections, the TQ-rclations were proven 
from the formula (JSJ. This formula contains the explicit factors w(z) and w(t), 
and one of them was incorporated into 77 (by defining i 7 JUjfc = w(t)Hj) and 
was associated to the nesting level. The other factor (w(z) = zS Xs) was de- 
composed into the characters of s-symmetric representations. 



16 Hirota relation can also be written in terms of the T-operators defined as T a>s (u) = 
T( a,s ) (u+a— s). In term of these T-operators, l|35jl takes the usual form T 0iS (u+l)T a , s (u—1) = 
T a +l «(^)T a _x s{u) + T a , s .|_i (u)T a s _i (it), more frequent in the literature. 

17 To prove the BR at arbitrary levels of nesting, one actually needs to rewrite l|86| l in the 
following slightly more general form, which is also equivalent to JH): 

det(^- fc W Iij («-2fc + 2)) jeJ 
W V W = 7 rr rjr-, 1 - fc - |J| , (36) 

det (*;-*) jeJ nlir wi(u-2*) 

i<fc<i./ 

which becomes the nested Bazhanov-Rcshctikhin formula after the limit Z{ — > l/xi is taken 
for i £ I . 
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But it is also possible to generate other identities from (JSJ, for instance by 
incorporating both w(t) and w(z) into II j, which gives rise to different nesting 
paths. Then, a careful rewritin g^! of ([5]) immediately gives the following well 
known QQ-relations between Q-operators: 

(xj - Xj) Qi(u-2)Qi, it j(u) = XiQi t j(u-2)Q Iti (u)-XjQij(u)Qi } i(u-2), (39) 

where i,j G I and i ^ j. These are relations among Q-operators in 4-cycles made 
of {Qj, Qi,i, Qi.j, Qi,i.j} in the Hasse diagram (cf. figure [T]). 

We can immediately solve these QQ-relations, which are actually a partic- 
ular case of the Pliicker identities for all the Q-operators to get the following 
determinant representation^^! 



det (x^ 1 k Qi.j(u — 2k]\ 
(nl J i7 1 Q/("-2fc))det (,•/ 1 ' ) 



^Au) = ,_,„„ . . N ■ (40) 

0<fc<|J|-l 



In particular, choosing 7 = gives the expression of any Q-operator in terms 
of K + 1 Q-operators, namely, the K single indexed Qi(u) operators describing 
the last level of nesting, and the u- independent operator Q$(u). More explicitly 
Qo(ii) can be defined by its action on the diagonal basis of the quantum space 

by 

K K , 

Q,(u)\e)^T s f°(u)\e) = 2 N l[n k \ ]J 1 - ^ |e), (41) 

k=l j=l, ^ k ^ 

where |e) = |e il:i2i ... jijv ) = \e lx ) (&\e i2 ) ® ■ ■ ■®\e lN _ 1 )®\e iN ) and n k is the number 
of j such that ij — k. 



4-5. Operatorial Bethe equations. In this subsection we derive the set of nested 
Bethe ansatz equations. 

From the QQ-relations (|39| one immediately sees that since Q/j(w) should 
be, by its definition, a polynomial of u, then (for i,j^I and i ^ j) 

Qi,i(u) | (xi - Xj)Q/(m - 2)Qj ) ij(u) + x i Q/ J (M)Q/ ji (u - 2), (42) 

Qi,i(u) | (xi - Xj)Q7(m)Q7,jj(u + 2) - a;iQ/ J -( u )Q/,i(' u + 2 ), (43) 

where P|P' denotes the fact that the polynomial P 1 contains the polynomial P 
as a factor. By adding XiQ/^u + 2) times the first line to XjQi^(u — 2) times 
the second line, one gets 

Qi,i( u ) I x iQi( u - 2)Qj.ij(w)Q/ /i (w + 2) + xjQi^Qi^j^ + 2)Q IA (u - 2). 

(44) 

18 In J5), let us put t = tj , z = ti, II = IIj^j Tfjyi wnere I C 3 and i,j £ I (i ^ j). Then note 
the following relations: w(z)w(t)II = w(z)II = TIjjj yj, w(t)II = Hj^jrn ■ 

19 The factor det (x'." 7 1 fe J jS j = Oi<j i — x j) m * nc denominator corre- 

0<*<|J|-1 

sponds to the denominator formula of the character of gl(\ J|). 
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This is written for the Q-operators, but when acting on particular eigenstates, 
the operatorial Bethe equations (|44|) become the usual polynomial Bethe equa- 
tions (gll) [28P29] (cf. eq. (68) of [22]) on the Bethe roots {u^} along a chosen 
nesting path. Indeed, now we know, by construction, that the Q-operators are 
polynomial^ of u and therefore for each eigenstate their eigenvalues are also 
polynomials of a degree Ki < N in u: 

Qi(u)=c I Y[(u-u[ I) ), (45) 
fe=i 

so that, substituting u = uj into an eigenvalue of (|44l) . we obtain the usual 
nested Bethe ansatz equations: 

for ke {1,2,..., 

(46) 

The Q-operators are formal polynomials Qi(u) = J2k=i uk whose coeffi- 
cients {cL } are also operators in the quantum space. On particular quantum 
states - the eigenvectors of the spin chain Hamiltonian, a part of them becomes 
zero, which explains the fact that the power of Q-functions - the eigenvalues of 
the Q-operators - can diminish on each step of the Backlund procedure. 

Given all solutions of the Bethe equations (j4"B"|) it is possible, in principle, to 
find all eigenvalues of the Q^. operators and then to reconstruct all T-operators 
using TQ-relations (see (|75j) ) together with the Hirota relation. 



Xj Qj(4 7 '° - 2)Q/,i(<' i) + 2)Q/,i,j(4 J ' i) ) 
*j Q I (4 / ' <) )Qj, i (4 /l<) - 2)0^(4^+2) 



5. Generalization to the supersymmetric case 

In the case of the gl(K\M) super-spin chain, T- and Q-operators are labelled by 
the 2 K+M subsets I of the full set 3 = {1, 2, . . . , K + M}. For any element of 3, 
we define the grading parameter: 

p b = for 1 < b < K, and p f = 1 for K + 1 < / < K + M. (47) 

Now the Q-operators are described by the colored Hasse diagram (see figure [5]). 
In the case of the gl(K\M) super-spin chain, the co-derivative can be defined 

by 

(48) 

where {e^} are generators of gl(K\M) in the fundamental representation and 
a matrix is expressed as <fi — ^2(/>j e ij7 an d g is a matrix in the fundamental 
representation of GL(K\M). 

20 This fact was missing in the analytic Bethe ansatz construction of [22| and appeared there 
only as a hypothetic ansatz for the solution of Hirota equation by the Backlund procedure. 
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Fig. 2: Hasse diagram for the Q-operators: g Z(2| 2) case. There are 2 4 = 16 Q- 
operators in the same way as gZ(4) case (figure fib]) . Qj and Q/\{fc} are connected 
by a solid line if k £ {1,2} (bosonic pk = 0), dashed line if k £ {3,4} (fermionic 
Pk = 1). 



As explained in [T], the properties of co-derivatives are exactly the same as 
in the bosonic case, including the expression of T-operators in terms of co- 
derivatives (at the zero th nesting level). The diagrammatics of the co-derivative 
is also the same as in the bosonic case, except the signs to be introduced into 
every permutation operator, to get V = *^2 a p{— l) Pa e^ Q ® e a p. In particular, the 
formula {8]) still holds in the super-case but ([7]) has to be substituted by 



w{z) = sdet — — = ^ l — — =Vz s Xs(g) = ^oo , 7vz ~ / aU r . 

l ~ z 9 Uf^ii 1 - * xj) ^ Ea= (-i) az a x { a >(g) 

(49) 

where (xi,--- , x K ,Vi, ■ • • ,Vm) = ,€k+m) are the eigenvalues of g £ 

GL(K\M) in the fundamental representation and sdet denotes the super-determi- 
nant. With slight generalizations of the definitions w.r.t. the bosonic case, all 
the supersymmetric TQ- and QQ-relations follow from ([HJ if we define T- and 
Q-operators in the following way: 



^7=n^-) (( " 



Bj=l[(l-Z j t j )(l-gt j )® N , 

jei 



(50) 
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(51) 



(52) 



wheren 5 = |Jn{l,2, ••• , AT}|, ancQnj = \ln{K+l,--- , K+M} |. In particular, 
the Q-operator for an empty set (j^Tj) now becomes 



K+M K+M / £ .n j -(-l) p i +p k 

Q ( M )|e)EETW( U )| e ) = 2 w [] {-l) Pknk n k \ J[ (l-f I 

fe=i j=i, ^ ^ fe - 



(53) 



In what follows, the indices i,j G {1,2,..., if} n I and fc, Z 6 {K + 1,K + 
2, . . . , K + M} (1 1 will correspond to the opposite gradings (we might call i 

and j "bosons" and k, I "fermions"), and we will use the notation I = 1 + K for 
I € {1, 2, ... , M}, so that & = yi. In (|5lT) . x\(gi) is defined through (|3"3"]) . where 
Xs(gi) is defined by the generating series 



10/(2) = sdet 



nr e i(i 

i>k 



z Vi) 



1 



9i 



3<K 



z Xj) 



s=0 



(54) 



That implies [50)17] that T< A} = if the Young diagram A = (A x , A 2 , A3 . . .) 
contains the highest weight \k+i > M (i.e. unless A is inside the "fat hook 

region" indicated in figure 10 in [I]), while for nested T-operators Tj A \ the 
same "fat-hook condition" holds, but the corner of the fat hook is displaced 
from (M,K) to {n f ,n b ), where n b = \I n {1, 2, • • • ,K}\, n f = \I n {K + 
1, • ■ ■ ,K + M}\. This means that the "hook region" decreases by one row or one 
column at each level of nesting. 

Then the TQ relations (|29]F1 become 

Tf (u)Q/j(u) = Tfj(u)Qj(«) - -r,T] /■„ + 2)Qj(u - 2) (55) 
TJ jf («)Q/(«) - TJ(«)Q Jtf («) - iftTJ-^u + 2)Q Jjf (« - 2) (56) 
for 1 < j < if and K + 1 < Z < K + M 



21 As before, |7 denotes Card(7) 

22 obtained from the master identity (0 by putting 77 = I lj^jr , t = tj for (|55 ft : 77 
77j,t = i-. for (SBJ. 
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The QQ-relation4E3 also become grading-dependent: 

(xi - Xj) Qi(u - 2)Q/.jj(u) = XiQi t j(u - 2)Q/.j(u) - XjQi : j(u)Qij(u - 2), 

(57) 

(xi - yi) Qj i(u - 2)Qi )f (u) = XiQ/(u - 2)Q 7 . t (u) - t//Q/(w)Q J i j(u - 2), 

(58) 

(Vi - Vm) Qj irni u ~ 2)Qj(u) = Z/zQj f(u - 2)Q/. I?l (u) - y m Q, I ? (tt)Qj )A (u - 2), 

(59) 

for i,j € {1,2,..., if} n I and f, m G {7f + 1, if + 2, . . . , K + M} R 7. 

These are relations among Q-operators in 4-cycles made of {Qj, Qi,i, Q/j, Qi,i,j} 
in the Hasse diagram fig. [U For example, in the figure eq. (l57|) corresponds to 
4-cycles made of 4-solid lines, ([55)1 corresponds to 4-cycles made of 2-solid lines 
and 2-dashed lines, (f59|) corresponds to 4-cycles made of 4-dashed lines. 

All these TQ- and QQ- relations are derived by choosing an appropriate 
II in ©, but, as explained in [3T], they could have been obtained from the 
bosonic relations by the so called "bosonization trick" : For instance, ()58|) can 
be rewritten as 

(x{ - yi) Qj(u - 2)Q J A; -(u) = XiQj^(u - 2)Q J)i (tt) - yiQ Js j(u)Qj(u - 2), 

(60) 

where J = I U {1}, which has exactly the same form as (f5T|) . up the the formal 
replacement J— > I, J\l — > 1,1. This is interpreted as the fact that adding a 
"boson" to the set I (which indexes Q or T-operators) is equivalent to removing 
a "fermion'Q- This trick can be viewed as a mnemonic rule, and in this con- 
struction of the T- and Q-operators, it comeJ^ from the power (— l) Pj in the 
definition (|5U|) of 77j. This trick means that the Hasse diagram for gl(K\M) can 
be "rotated" by putting Qk+i,K+2,— ,k+m on the top row and Qi,2,-- ,K on the 
bottom row, so that the QQ-relations take the same form for all the facets of 
the modified diagram. 

In the same way as for the gl(K) case, the TQ-relations (|55ll56"|) can be written 
as the generating series expansion ((75)) (see [2T)1[2"2"] at the level of eigenvalues). 
We can also generalize the Wronskian expressions (j4H]) for Q-operators to all 
the T-operators at all levels of nesting. They do not differ in the form from the 
relations found in [2"5] . 

The Bethe ansatz equations for gl(K\M) [5SJ are obtained in the same way 
as for the gl(K) case: there are two "bosonic" Bethe Ansatz equations (BAEs) 

23 obtained from the master identity (0 by putting 77 = 77-j— ^fjj, t = tj,z = ti for H57jl : 
77 = n w{i} , t = tp z = ti for (58}; 77 = 77j, t = t A , z = tj for |g§). 

24 As explained in |22| . the same linear system describes the addition of a column to the "fat 
hook" or the removal of a line (see (41) and (42) in 20 22] ). In our notation, the former case 
is the addition of a "fermion" to the set I, while the latter is the removal of a "boson" . 

25 In (150152ft we can see that adding a "boson" to the set 7 multiplies 77j by w(tj), while 
adding a "fermion" divides it by w(tj). As the "Master Identity" is only sensible to the addi- 
tion/removal of w(tj) factors, it is not surprising that the QQ-relation obtained from "master 
identity" satisfies this "bosonization trick." 
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and two "fermionic" onesB The first "bosonic" BAE follows from (|57p and is 
unchanged with respect to the section 14.51 : eq. (|44| on the level of the operator 
and eq. (l4"6")l on the level of the eigenvalue. The second bosonic BAE is obtained 
by isolating Qj ? in (jSUj) which gives 

Q 7 i(u) | ymQj irh {u - 2)Q J (u)Q J i(u + 2) + y ; Q / f A («)Q/(« + 2)Q, r (u - 2), 

(61) 

which is, at the level of eigenvalues, the equation (|8ip . equivalent to the equation 
(69) of [22] (up to the permutation on the indices (the Weyl group symmetry)). 

On the other hand, the "fermionic" BAEs can be immediately obtained from 
(|58p in the form 

Qi,i( u ) I XiQi(u - 2)Q I i i (u) - yiQj^Qj ^u - 2), (62) 

Qi,f(tt) | x l Q I {u)Q I iJ (u + 2) - yiQi(u + 2)Q I iJ (u). (63) 

In terms of eigenvalues, the Bethe equation (]62l) (resp (J63J) ) is written as (|82|l 
(resp ([55)1 ) in the appendix [D] 

6. Conclusions 

The co-derivative formalism and the master identity ([HJ , together with the defi- 
nitions ([22]) . (|32 |) of nested T- and Q-operators proposed in this paper can serve 
as an alternative approach to the quantum integrability, rather different from the 
popular algebraic Bethe ansatz (see for example [32 and the references therein). 
It allows to complete the whole procedure of diagonalization of transfer-matrix 
of the inhomogeneous twisted gl(K\M) (super)spin chain, all the way from its 
construction from R-matrices obeying the Yang-Baxter relations and till the 
nested system of Bethe ansatz equations, directly in terms of the operators act- 
ing on the quantum space. The master identity © presented at the beginning 
of the paper and generalizing a similar identity from 1 is the basis of this ap- 
proach, encoding all possible operatorial QQ- and Baxter's TQ-relations at every 
step of nesting, or of the operatorial Backlund flow, generalizing the operatorial 
Backlund transformations of 21,2G\. Remarkably, the master identity takes a 
bilinear form with respect to the gl(K\M) characters, or their generating func- 
tions. Since the characters can be viewed as the tau-functions of KdV hierarchy 
(which is of course a reduction of KP hierarchy) one can speculate that this 
identity is simply a particular case of the general Hirota identity for the KP 
tau-functions, with r„ = itr g n playing the role of the KP "times" . It would be 
an interesting relation between the quantum and classical integrability, showing 
that, paradoxically, the former is a particular case of the latter. 

It would be also interesting to generalize our approach to the case of non- 
compact representations of gl(K\M) in the auxiliary space, following the ob- 
servations made in [53113 lll3"4"] for the characters and Q-operators for {7(2, 2 14). 
This might teach us how to deal with one of the most interesting integrable 

26 We will call "bosonic" (resp "fermionic") the Bethe Ansatz Equations having two free 
indices of same grading (resp of opposite grading). For instance, in Ij61|l . the two free indices I 
and m have the same grading, hence this Bethe equation is called "bosonic" . 
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physical systems, J\f = 4 SYM theory and its AdS dual - the Metsaev-Tseytlin 
sigma-model having the PSU (2, 2|4) global symmetry. In general, the Y-systems 
for sigma-models and their Wronskian solutions 8,9 might be also an interest- 
ing subject for their operatorial generalization in the quantum (physical) space 
and might give us an interesting tool for the study of the spectrum of excited 
states and shed some light on the formulas for the energy of an excited state 
conjectured in the literature for relativistic sigma models [35 l l36 [ l9] and for the 
AdS/CFT [37]. It would be also interesting to generalize our master identity to 
other symmetries of the spin chains, where the Bazhanov-Reshetikhin-type re- 
lations are also known, to other than fundamental irreps in the quantum space, 
as well as to the trigonometric and elliptic R-matrices. 
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Note added 

When finishing this paper we learned about the results of [H] which deal with 
the same problem, the operatorial formulation of the Q-operators and of the 
TQ-relations. The objects studied in that forthcoming paper are the same, but 
the formalism is radically different from ours. 



Notes for version 2 

We corrected in this version [arXiv:1010.4022v2 [math-ph]] many minor mis- 
prints of the first version, and added some details and explanations. It is pre- 
cised that 11(g) in the master identity (jHJ) has a form of det /(<?), where f(z) is 
an arbitrary fixed function, rather than an arbitrary class function of the twist 
matrix g. We also removed a false TT-relation (eq. (4.27) in the version 1). 

Most importantly, we give in this new version a concise proof of our master 
identity - the basis of our approach. 



Notes for version 3 

Minor misprints are corrected. 
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A. Diagrammatics of co-derivatives 

A.l. Co-Derivatives and characters. The action of the co-derivative on charac- 
ters and their generating function is explained in pQ. For instance, we can write 

D ® D ^ x) = _£E_®_e_+^ ('^_ g) = HI + V) «,(*), 

1 — gx 1 — gx \ 1 — gx 1 — gx J \ 1 1 •' k / 

(64) 

where 1 (resp |) stands for j^- (resp t^— ), so that I I = j^- ® rj^- and 

XI it Q-E 9 X 11 9 X 9-E 

= ~Pi 2( 1 \~ ® 1 jj"L ) , where "Pi 2 denotes the usual permutation between the 

1 st and 2 nd quantum spacJ^l. The general expression for D® N w(x) is given by 
the formula (4.11) of [1. (formula (30) in the arXiv version), and it represents 
the sum of diagrams corresponding to all possible permutations, with dashing 
all the lines going from lower to upper nodes and directed to the right. 

In terms of characters of the symmetric tensor irreps Xs(<?)) the equation (|64l) 
reads 

D®L>Xt=(J\ + )Qxt where J (resp I) stands for ^ g *lo t (resp Y _^ e -a t )■ If 

we identify x = e~ dt there is no ambiguity between these two definitions of the 
diagrammatics. 

This diagrammatics can also be extended for the inclusion of the parameters 
Ui, i.e. 

( Ul +D)®(u 2 +D) w(x) = (II + II + H + II + XH^' ( 65 ) 

where | stands for u\ or u 2 , according to its position (for instance, | J = ml ® 
and If = ^®lu 2 ). 



A. 2. Co-derivatives of products. In this paper, we often use the co-derivatives 
acting on products, defining the quantities like D® N 77 ■ w(x) for an arbitrary 
class function 77. Then each co-derivative can act either on 77 or on w(x), and 
at two spins, for instance, the Leibniz rule gives 



D(g>D n w{x) 



D 



d n 



d®d n 
b®i n 



w{x) 
w{x) + 



n 



b w{x) 



i®d n 



n 



n 

'jj(x) 



I® 73 to (a:) 

+ 



n 



(66) 

73® I w[x) 
D®D w(x)] (67) 



11 

'Jl}(x) 



n 



n 



n 



n 

v(x) 



(68) 



n 



(69) 



27 For the supergroups, 'Pi, 2 is replaced by the super permutation which differs from the 
usual permutation only by certain signs (see pQ). 
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The equality between ()66|) and (|S7|) is just the Leibniz rule, while defines 
the graphical representation of each term of (|67|) . Each black dot stands for a co- 
derivative, acting on what lies on its right on the same horizontal line (horizontal 
lines actually hide auxiliary spaces whose characters contribute to 11(g) [resp 
w(x)] ). The vertical lines correspond to the quantum space on which the whole 
operator acts, and the crossings without dots stand for I. When operators are 
multiplied in the quantum spaces, they are represented one above another. For 

- n 

- w(x) 



instance 



• D n 



t)®I w(x) 



The last expression gives a shorter representation of (|S7|) . where an 
implicit multiplication is taken between the blocks ^ +^ ■ The generalization 



to an arbitrary number of spins N is straightforward - one has iV such blocks 
instead of two. 

One can also see that due to the relation Ddet(g) = Idct(g) and to the 
Leibniz rule, we have very generally 



)(u 1 + 2D) ndet{g) a 



= det(g) a 



){u! + 2a + 2D) 77 



A. 3. Commutativity of all T- and Q-operators. Everywhere through this paper, 



<8>£i(«i 



D) 77(g) , where 11(g) 



the quantities of interest are of the forrrp 

is a class-invariant function of g (a symmetric function of its eigenvalues). 

In [TJ, the particular case when 77 is the character of a representation {A} 
was studied, and in particular it was proven that the T-operators commute for 
different representations {A} and {//}, yielding the relation 



N 



(ui + D) x{\}(g) 



N 



)(m+v + D) X{p,}{g) 



1=1 



0, (70) 



where \A, B\ = A ■ B — B ■ A, and Uj, v S C. Then, by writing an arbitrary class 
function as a linear combination of characters, one gets 



jv 



)(ui+D) 77 



N 



)( Ul + v + D) 77' 



= 0, 



(71) 



which holds when 77, 77' are two arbitrary class functions of g. 



B. Degree in u of the polynomials Tj(u) 

We already claimed that the l.h.s. of (TT2j) only has simple poles. Now that we have 
described diagrammatic rules for the action of co-derivative on w(z) functions, 
we can make this statement more explicit by writing the matrix coefficients of 
the Q-operator: 

28 Most often, the operators of interest were actually of the form ^'^Lj^f i + 2 D) i7(g)j , 
which is equivalent, after the rescaling m — > 2 U{ and 77 — » II/2 N , to ®£Li ( u * + ^0 ^(s) • 
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These matrix elements are given by a generalization of the formula (4.16) 
[(35) in the arXiv version] of p], which reads (see also (|65|)f^l 



){2 + Ui + 2D) w{t) 



N 



En 



2(ff t) 



6(a(i)-i~l) \ fe ' 



w(t), (72) 



where the sum is taken over the permutation group Sn of order N. At this 
point, we can already notice that the only non-zero terms in the sum are the 
permutations such that Vi, fcj = l a (i), since we are working in the diagonal basis. 
The same property holds for T-operators, and it implies for instance that the 
number of spins pointing in each direction \ek) is a conserved quantity. 

After multiplication by (1 — gt)® N , the left-hand-side of (fT2|) has only simple 
poles from w(z)\ z ^ 1/xk : 



[i-gtf 



)(u l + 2D + 2) w(t) 



fel ,/E2 ! "'^J v 



N 

E II 5 t m ( u ^°d) (1 - **« *) + 2 fe t) e ^-*-V) w(t). (73) 



In the definition of Qj(u), we take the limit £ — > 1/xj, and we see that 
Qj('«)|efc l! fc 2i ...fc JV ) is independent of all u n such that fc„ = j. As a consequence, 
the degree in u of each eigenvalue of Qj(u) is equal to the number of spins 
pointing in the directions \ek}k^j in the corresponding eigcnstatc. 

Moreover, the generalization to an arbitrary T or Q-operator is very simple: 

T| A ^(M)|e; l! ; 2 .... t i N ) is independent of Ui for alH such that U G /, and the degree 

in u of each eigenvalue of t| A ^(m) is equal to the number of spins pointing in 
the directions \ej)j^i in the corresponding eigenstate. 



C. Generating series of symmetric T-operators 

The TQ-relation (|29p can be re-written as an expression for the generating se- 
rieS 2H/ (u, z) = J2T=o zST i ( u ) of symmetric T-operators: [2ll[20|[22|[7] 

W*^!^ 1 -'^ 1 *)**' 4 ,74) 

where we used a shift operator e 29 "/(u) = f(u + 2) for any function f(u). Let 
us then fix a chain of subsets of the full set: 3 = Ik D Ik-x D • " D Iq = 9, 



Here, the function 8(n) is equal to 1 [resp 0] if n > [resp n < —1]. 

From Xs{g<l)) = <5s,0j we can see that W${u, z) = Q0(m), which we will use to get H75I I. 
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where I k = {ji,j 2 , • • ■ ,jk}, h \ h-i = {jk}, k = 1, 2, . . . , K. Then from (H), 
we immediately get that 

2H/ fc («, «) = Ok • Ofc-i • • ■ Oi ■ Q0(«), (75) 
where 0fc = fl-^%i^^V X 



n=0 



which expresses the T-operators for symmetric tensor representations in terms 
of the Q-operators Q/^. In particular, taking the coefficient of z in (|75|). one gets 
the usual relation for the fundamental representation s = 1: 

K (U) = Q/ * (U) £ -q^) Q/m _ i(u) ■ (76) 

Note that this T-operator ([?!)]) has the same form as the T-function in [25] 
obtained by the nested Bethe ansatz, if the Q-operators are replaced by their 
eigenvalues (Q-functions). 

The same generating series can also be written in the gl(K\M) case for a fixed 
chain of subsets of the full set: 3 = Ik+m 3 Ik+m-i D • • • D Io — ®i where 
h = {ji, h, ■ ■ ■ ,3k}, h \ 4-1 = {jk}, k = 1, 2, . . . , K + M: d still holds for 
j < K, while for I > K + 1 gives : 

ttW"'*' - { l - "'^r^'j (77) 

so that (JTSJ) becomes for gl(K\M) : 



W Ik («, z) = O fc • fe _i ■ • ■ Ox • Q (tt), (78) 



Q/ fc -x("-2) 2a \ 1 Q/» 



where O k = ( 1 - x Jfc ^~ ze 2d " ) ^ if jfc < K, 



which gives for instance the following generalization of (|76[) : 

T i r s n , , y^, nPifc , Q/ fc ->-2(-l)^)Q 7 > + 2(-l)^) 

(79) 

Note that the eigenvalue of (f79|) coincides with a traditional form of the T- 
function from the Bethe ansatz [35]. There are (K + M)\ ways to chose the 
chain {ife}S) M , ^" u ^ f79"l) does not depend on this choice. This is the (super) 
Weyl group symmetry of the T-operators. 
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D. Bethe equations at the level of eigenvalues 

As explained in section l4~Sl the operator equation (|44p can be written as a Bethe 
equation (|4l)| on the roots of the eigenvalues Qi(u) of the polynomial operators 
Ql(u) for any eigenstate. 

The same can be easily done in the supersymmetric case so that ([44)) . 
([621) and (|63|) imply that for any eigenstate, the eigenvalues of the Q-operators 
satisfy the conditional : 



^Q / (4 / ' l) )Q/,(4 / ' l) -2)Q7«(4 7 ' i) +2) 



1 in Q^(4 J ' ) +2)Q 7 >P-2)Q UA (4^) ral , 

— 1 = ; ^ — — ; lor 1 < k < K T f , 81) 

^QH4^)Qz,r(4 / '° + 2)Q Z;U (4^-2) 

Q/(4 f,0 -2)Qw(»i f,<) ) , , <h<K rs9 , 

1= «QK^)Q^-2) for l<„^,( 82 ) 

1 = E QK^^2)Q JAf (u^) forl<fe< ^, (83 ) 

^q,(4^)Qw(4^+2) " " u 

where (4 j ■ ■ ■ , vSS ') are the roots of the polynomial Qi(u). 

These equations (|80M83|) arc equivalent to the traditional form of the Bethe 
equations [23], respectively, [cf. eqs. (68), (69), (71) and (70) in [55].] 



E. Proof of the master identity ((SJ 

In this appendix, we will give a proof of the master identity (JSJ based on the 
so-called Bazhanov-Reshetikhin formulgP^I. 

Consider the main object of our Master identity the generating operator of 
the transfer matrices 

W ; (u) = + 2D) ]J w(z k ), (84) 

i fc£j 

where zjeC and I is any subset of Z>o- We assume Zi ^ Zj for any i,j 6 Z>o 
such that i ^ j. Then the master identity ([SJ can be rewritten in a form of a 
Pliicker identity (or Jacobi identity): 

( Zi - Zj)Wi, itj (u + 2)Wj(u) - ZiW Itj (u)W Iti (u + 2) - zjWjjfo + 2)Wj,,(tt), 

(85) 



31 Like in the section[5] we use here the convention i,j € {1,2,..., K}, and denote by I = l+K 
and m = m + K for 2, fe 6 {1, 2, . . . , M} for some indices of fermionic grading. 

32 We thank Anton Zabrodin who proposed us the idea to use the Bazhanov-Reshetikhin 
formula for the proof of our master identity. 
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where i,j G Z>o, i, J ' ^ I, i ^ j. It can be solved recursively in the same way as 
the QQ-relations (|40jl) giving 

det (z"- fe W,(w-2fc + 2)) 
W/(«) = / , . ;yi <^ fc <" , (86) 

det (-r fe ) 1 <„<„n^ 1 ^-2fc) 

Here (f>(u) = Wj(u) = rX/_i u i = Q^C")) an d we consider (without losing 
generality) the case I — {1, 2, . . . , n} for any finite n £ Z>o- 

Hence, to prove the master identity ([£} all we need is to prove this determinant 
formula (j8"rJj) . Let us expanc0 the left-hand side of (|5oT) multiplied by a factor 
with respect to {zk}'- 



n 

Wj( u ) det («?-*) TT = det (zr k ) w du) 



i=l 

oo oo oo n 



= <^,(u i+ 2D) £ sgn(<7) 2 E ••• E II W^r"^' (87) 

i <7£S n mi= 0ni2— ni n — Ok— 1 

where the surrF^I is taken over the permutation group on the set {1, 2, . . . , n} 
and sgn(er) is a signature of the permutation a. The coefficient of Ofc=i z k k ^ or 
any set of integers {Xk}kei in |57|) . due to the Jacobi-Trudi formula is 
nothing but the transfer matrix T\ from ([6]) 

n 

+ 213) £ sgn(a) [] XA fc+CT (fc)-fc(5) = 

= (g)(« i + 2 J D) det (xA 3+fc -,(.g)). (88) 

i 

Now let us expand the right hand side of ((SB]) (times the same {z/cj-dependent 
factor): 



i<ss eB (^"* w ^ tt - 2 * +2 ))n^" B 

OO OO OO 71 

= £ sg n( CT ) £ e ••• E n Tmfc ( M - 2 ^)+2)zr~ aW+fc - ( 8Q ) 

crGSn mi —0 rn2—0 m n — Ok— 1 

The coefficient of flLi 4* in ® is 



V sgn(o-) TTT Afe+<T W- fc (u-2CT(fc)+2) = det (T^ +k ~ 1 (u - 2k + 2)) . 

(90) 



1<] 

cr€5„ fe=l 



33 Eq. I I87C is an operator analogue of a generating function of the T-functions for any Young 
diagrams in eq. (2.44) in |23] . 

34 The sum over {m^^^j can be taken over any integers since Xmig) = if m < 0. 
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Therefore, the proof of (fSijj) is reduced to the following identity 
*r*, ^ , > \ deti<, fc <„ (T x i+ k -Uu-2k + 2)) , 

®(„, + 2D, A (x^-M) = ' n L;, ( J M) ^ <»D 



where the l.h.s. is precisely the T-operator T^ x ^(u). We recognize here the 
Bazhanov-Reshetikhin formula proven in 1,. This proves the formulas ((55|) - ((55|) . 
and hence the master identity (f5|)P^I 



F. Co-derivative and "removal" of eigenvalues 



The co-derivative of Xs(gi) a priori does not have such a simple expression (as 
(|64l) in terms of diagrams) as the co-derivative of Xs(ff)- We will see in this 
subsection how to define the action of co-derivative on Xs{gi), and then we will 
see how to compute the corresponding T-operators. In particular we will prove 
the relation (|2"g)l . 



Definition. 

The action of co-derivatives on wi(z) (introduced in (|23jl ). can be defined by 
means of equation ([5]), provided we specify what Xj is at the point e^' e g. The 
most natural definition is based on the fact that Xj is the j eigenvalue of g, or 
in other words the j th root of its characteristic polynomial. In this sense, Xj is a 
function of the group element g: Xj — Xj(g). In particular, Xj(figf2^ 1 ) — Xj(g) 
for any similarity transformation. 

If g is a diagonal matrix, it is immediate to see that the contribution of 
the non-diagonal-elements of the matrix e^' e g to the characteristic polynomial 
det (AI — e^' e g) is at least quadratic in </>. This means that at the point e^' e g, 

Xj is equal to (e*' e <?)^. to the first order in (f>. As a consequence, we get D^Xj = 

jiQj ~ ^j^ 3 ji x oi so tnat Dxj — PjXj, where the projector to the eigenspace for 
the j'-th eigenvalue Xj is Pj — ejj in this case. 

More generally, if g = Q~ l qQ where g is diagonal and Q is an arbitrary 
similarity transformation, then 36 ! we obtain 

(92) 



This exactly means that for a non-diagonal matrix <?, D Xj = PjXj, where 
the projector to the eigenspace for Xj has the form Pj = Q~ l e.jjfl. 



35 Note that both sides of (191ft are antisymmetric w.r.t. the set (Ai — 1, A2 — 2, . . . , A n — n). 
Hence, when the determinant is non-zero, we can always relabel A^'s in such a way that the 
highest weight components satisfy the usual inequalities: Ai > A2 > • • ■ > A n > 0. 

36 In (192 ft . we mainly use the relation (2.6) [equation (12) in the arXiv version] of [l]. 
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Computation of T- operators. 

The claim which was already given in (|28|) is that the computation of T- 



operators is done by commuting a factor 



u T (z) 



to the left of the co-derivatives. 



We will show that in the definition (l?7|) (or equivalently (JUJ)), the multipli- 
cation by Bj introduced in (f^Uj) allows to commute any function of Xj (where 
j G /) across the co-derivatives. As a consequence, the right hand sides of (|27l) 
and (j2"51) arc equal, giving 



5> s t?( U )= [Ha 



lim B 7 - 



JV 



)(u i + 2D + 2|J|) w(z)7T 7 



(93) 



where the r.h.s. can be easily computed by diagrammatic methods. 
In the case of 1 spin, this is checked by computinj^l 



lim (1 - gt) (u + 2D),Xjl 



(94) 



We see that the key point in this commutation is the multiplication by lini* , j_ (1 — 

gt) = (1 — g/xj), which cancels the terms in D Xj — XjPj due to the property 
(1-9/^)^=0. 

At N = m + n spins (m G Z>o, n G Z>i), the analogous relation is 



where C mj „ = (1 - .9/%) 



(95) 



(m + 2D) 



\i— m+2 



and it is proven by the recurrence over m. For m = 0, this follows from (|94j) . 
Let's show how C m +i, n cancels under the assumption that C m ^ n = for all 
g G GL(K) and any itj G C, j = 1, 2, . . . , m + n. Then for any uq G C, one can 
calculate: 



= ((1 - g/x 3 ) ® . (( Uo + 2 £) ® C m ,„) 

=C; i+ i, i + 2((l- 5 /x J )®I 8(m+ " 



(96) 



Here, [A, B] denotes the commutator AB — BA. 
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wheri 



38 



^m+Ln — (1 9/ x j) C 



Bi 



m+l,n 



(rn+n+1) ^/ 



m+l,n> 



)(«i + 2£>) 



(Urn + 1 + 2D),XjI 



m-\-n 



(u< + 225) 



y? — m + 2 



(97) 



and due tco D ® g/ x j = V ■ (1 ® g/xj) — Pj ® <?/:£;, the second term in ([96 
(multiplied by 1/2) can be expanded to get 



n' - 



: ((i - g/xj) 
(p.-g/xj)t 



(m+n) 



(m+n) 



f m+n 



' m+n 



\ k=l 
/ m+n 



'(l-Sf-g/xj) ]-{I®B m>n ) 



'm+n 



1(1-^-5/^-) •(I®B m ,„) 



. k=l 



and the first term becomes 



E^ofc ■ (l® fc ® (ff/xfc) ® • (l ® (1 - fl /a; fc )® (m+n) ) • (I ® B m , n ) . 



whicrjjfl is zero because it contains C m ,n- The second term is also zero because 
it contains (1 — g/xj)Pj. This completes the proof of the fact that C' m+l n = 0, 
from which C m +\,n = follows. 

As a consequence, we can indeed commute the factors gjr^ to the left of all 
co-derivatives in ([27)1 and get the relation 



38 In @6j, we used the Leibniz rule lW(l - g/xj)® m+n • B m ,n) = [i? <g> (1 - g/xj) 9m+n ] ■ 
(I ® B m>n ) + (l ® (1 - ff/a:j)® m+n ) ■ [D ® B m ,„] . 

39 Here, the projector Pj on the the j th should not be confused with the permutation operator 



V between the quantum spaces. 

40 Here, we use the fact that ((1 - -jL) ® E®( m + n )) -V , k = V 0>k -I®*® (1-fl/a 



n+m — fc 
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